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Abstract 

We consider cosmological models with dynamical dark energy (dDE) coupled to 
cold dark matter (CDM), while simultaneously allowing neutrinos to be massive. 
Using a MCMC approach, we compare these models with a wide range of cos- 
mological data sets. We find a strong correlation between this coupling strength 
and the neutrino mass. This correlation persists when BAO data are included 
in the analysis. We add then priors on v mass from particle experiments. The 
claimed detection of v mass from the Heidelberg-Moscow neutrinoless double-/? 
decay experiment would imply a 7-8cr detection of CDM-DE coupling. Simi- 
larly, the detection of v mass from coming KATRIN tritium /? decay experiment 
will imply a safe detection of a coupling in the dark sector. Previous attempts 
to accommodate cosmic phenomenology with such possible v mass data made 
recourse to a w < —1 eoS. We compare such an option with the coupling option 
and find that the latter allows a drastic improvement. 



1. Introduction 

Today we have a standard model of cosmology, the so-called ACDM model, 
providing an excellent fit to all cosmological data. This model tells us that dark 
energy (DE) and cold dark matter (CDM) account of ^ 75 % and 20 % of the 
present cosmic energy budget, respectively 'il]. It is then quite embarrassing 
that we fail to understand the nature of both DE and CDM. Furthermore, if 
DE is a smooth cosmological fiuid with an equation of state w = — 1, our model 
is troubled by two fundamental questions related to its relative and absolute 
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density, the coincidence and fine tuning problems: Why did DE start to domi- 
nate just when structures had time to form? Why is its density ^120 orders of 
magnitude smaller than the (naively) expected quantum vacuum density? 

One way to circumvent the coincidence problem amounts to introducing a 
coupling between CDM and dynamical DE (dDE) [3, |3] ■ The energy transfer 
from CDM then allows DE to comprise a significant fraction of the cosmic energy 
budget over a large part of the cosmic history. Unfortunately, switching on a 
coupling apparently worsens the fit of cosmological data. 

More recently, however, it has been noted that neutrino masses (ttIi/) and 
a CDM-dDE coupling affect cosmological observables in opposite ways [1, Q . 
However, not only our ignorance about m^, softens the constraints on the cou- 
pling, but a much more puzzling effect arises: models with ^ and coupling 
appear (slightly) favored in respect to m,^ ~ uncoupled models. 

Neutrinos are abundant in the Universe, second only to photons when it 
comes to number density; v oscillation experiments tell us that they are massive 
and measure the mass split between i/-mass eigenstates, so that the largest 
split O.OSeV from atmospheric ly's) is a lower limit on the heaviest v mass. 

Particle experiments have placed various upper limits on the absolute v 
mass scale. The Mainz and Troitsk experiments, measuring the end-point of 
the electron energy distribution in tritium /? decay, gave a 95% C.L. limit < 
2.0eV 0. A further controversial detection of an absolute v mass came from 
the Heidelberg- Moscow (HM) experiment; on the basis of its outputs, a part 
of the experiment team claims a 3cr lower limit m/j^ = (0.2 — 0.6)eV 0, [13] 
(the mass measured by neutrinoless double /3 decay (Oi^/3/3) experiments, m/3/3, 
is however a different combination of mass eigenvalues than the mass mp, from 
tritium /3 decay). 

In 2011, the experiment KATRIN [11], also studying tritium (3 decay, is 
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expected to start taking data. With a prospected sensitivity of cr^2 « 0.025eV 
it should be able to detect rui, values in the range of the mpfj claim. 

At present, the best upper limits on the v mass scale come from cosmology, 
yielding Em^ = < 0.2eV and M, <1.5eV (at 95% C.L.), depending 

on data sets and cosmological models 12 , 13 , 3, IE 3 H: S 

All this is true if the possibility of a dDE-CDM coupling is neglected. 

In this paper, similarly to what is done in Q, we will allow for such a 
coupling, reporting for the first time results when baryonic acoustic oscillation 
(BAO) data are taken into account. Then we consider the v mass limits set by 
a part of the HM collaboration (KKDC claim, herafter), and impose this as a 
prior on M^. Finally we shall investigate how a mp detection from KATRIN 
would affect limits on coupling. 

The possibility that forthcoming neutrino experiments yield mass values in 
apparent conflict with cosmic data had been considered by various authors. The 
most promising option, perhaps, had been discussed by [18|, , who found that 
allowing w, the state parameter of DE, to take values < —1 eased such constrast. 
Here we compare this option with the coupling option and find that the latter 
one leads to a drastically better agreement between terrestrial and cosmological 
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measures. 

In the following section we discuss our cosmological model and outline the 
observable effects of the CDM-dDE coupling and massive i^'s. In Section [3] we 
discuss the experimental bounds on the absolute v mass scale. The data and 
methods used are presented in Section SI while Section [S] is devoted to reporting 
and discussing our results. In SectionlHlwe summarize our findings and conclude. 

2. Cosmological model 

Our models differ from the standard ACDM in three different aspects: (i) 
DE is a self-interacting scalar field rather than a cosmological constant A. 
(ii) A linear dDE-CDM coupling is allowed, (iii) We allow z^'s to be massive. 

We shall consider the Ratra-Peebles (RP) potential [20[ and a SUGRA self- 
interaction potentials |2l| , reading 

V{4>) = A'^+Vr , Vicf,) = (A"+Vr ) exp(47r</.V™p) , (1) 

respectively; they allow tracker solutions for any a > 0. For both potentials, 
once a and A are assigned, the DE density parameter VIo^de is uniquely defined. 
In our fitting procedure, however, we use ^o,de and A = log(A/GeV) as free 
parameters. 

Limits on these models without coupling between DE and CDM have been 
studied in For most cosmological parameters, WMAP5 results lead to a 

slight narrowing of error bars, in respect to WMAP3. In the case of A, however, 
[22] find a significant shift downward of the 2-a upper limit on the energy scale 
A. In the SUGRA case, in particular, only A < —3.5 is allowed. Such small 
values are well below the range motivated by particle physics. Therefore the 
physical appeal of the SUGRA potential is spoiled. 

Following the procedure in refs. 0, we assume a linear coupling between 
the DE and CDM energy components, which leads to the coupled equations 

4> + 2-^ + a^V'{(f>) :^Ca^Pc (2) 
a 

Pc + 3-pc ~C(j)pc ; (3) 
a 

here a is the scale factor, pc is CDM density, C denotes the coupling strength. 
We then define the dimensionless coupling parameter 

/3 EE \/3/I67r nipC, (4) 

used as a free parameter for our cosmological models. For a more thorough 
discussion of the effects of coupling on dDE and CDM evolutions, see, e.g., [5|. 

Let us however outline that, when the [3 degree of freedom is opened, A 
values as large as 30GeV become allowed, at the 1-a level, while, at the 2-a 
level, no significant constraint on the energy scale A remains. Even for the RP 
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potential, for which a limit A < —8.5 held, in the absence of coupling, values 
A ~ — 2 become allowed. 

The effects of massive neutrinos in cosmology have been studied thoroughly 
for many years, and we refer to [31 for an extensive review of the topic. Cos- 
mological observations are mostly sensitive to the sum of i/-masses, Mj^, related 
to the V density parameter by the relation fl^/i^ = My/93. 14eV. 

If My < 4.5 eV, most i/'s were still relativistic at matter-radiation equality 
and act then as radiation, so postponing the equality compared to a model where 
all the DM is cold. When keeping the total DM density constant, increasing 
rijy shifts the peaks of the anisotropy spectrum Cg to lower £ and boosts their 
heights. 

When it comes to the matter power spectrum P{k), the main effect is ly free- 
streaming from small scale fluctuations, damping P{k) for large wavenumbers k. 

The combined effects on Ce and P{k) spectra, when c omp ared with obser- 
vations, lead then to the limits outlined above [H, [13, [H, 13, 3 Q These 
effects, however, as outlined in @,@], are almost opposite to those of dDE-CDM 
coupling, for both Ci and P{k). A Fisher-matrix analysis then shows a strong 
degeneracy between the /3 and M^, parameters. Pinning down one of the pa- 
rameters by some other means clearly results in improved limits on the other 
parameter. 



3. Neutrino mass bounds from earth based experiments 

Two different ways are being followed to measure the absolute scale of v 
masses: Qvj3l3 and tritium (3 decay experiments. 

The Qv^P process can only occur if z/'s are massive Majorana spinors, i.e. 
they coincide with their own antiparticles. From the measurement oiT^J^ (Oi'/?/? 
half life) one can then deduce an effective mass m^^ = U^^rrii, being 

2 ''^e /r\ 

= r — to^- 

Here ?7ei is the PMNS v mixing matrix, mi are mass eigenvalues, me is the elec- 
tron mass. Cmrm the nuclear matrix element relevant for the nuclide considered, 
is the main problem with Ovfi/i, because of its large theoretical uncertainty. 

The most sensitive limits, up to now, were derived from detectors enriched 
in ^"^Ge. Using such nuchde, the Heidelberg-Moscow (HM) [H and the IGEX 
Hi gave the hmits T^^^ > 1.9 x lO^^y and T^J^ > 1.6 x lO^^y, respectively. 



However, a part of the HM collaboration published results claiming a more than 
5(7 detection of a signal, giving an effective neutrino mass of m^p = (0.2 — 0.6)eV 
(3cr hmits) [i,[l3 (which we refer to as the KKDC claim). 

In view of nuclear matrix uncertainties, this claim is not in contradiction 
with another experiment, CUORICINO 25|, based on ^'^'^Tc, placing an 

upper hmit of m,3^ < {0.19 — 0.68}eV. 

Both ^■''^Te and ^^Ge double-/3 decays are however still being investigated 
through the CUORE and the GERDA experiments, respectively. The CUORE 
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experiment [2^ is running in the Laboratori Nazionali del Gran Sasso (LNGS), 
where the GERDA experiment 27[ is also being placed. 

Tritium /3 decay is a different road to the absolute neutrino mass scale, where 
one accurately measures the energy distribution of the outgoing electron, and 
uses this to infer an effective electron neutrino mass, to| = J^i \ Uei \ rnf. 

Currently, the best limits on from tritium /3 decay come from the Mainz 
and Troitsk experiments, yielding an upper limit of < 2.0 eV (95% C. L.). 
The much more sensitive KATRIN experiment [lH is scheduled to start taking 
data this year. When completed, KATRIN is expected to reach a sensitivity of 
a^2 « 0.025 eV^, and thus be able to confirm the KKDC claim. 



4. Data and methods 

We use a modified version of the public program CAMB [28^ to calculate 
CMB and matter power spectra. A modified version of the public MCMC engi- 
ne GosmoMG [1^ is used to generate confidence limits on the parameters: {wfc, 
ujc, T, Us, In 10^^ As, logipA, P, M^}. Here uJb,uJc are the reduced density 
parameters of baryons, GDM; 9 is the ratio of the sound horizon to the an- 
gular diameter distance, r is the optical depth, Ug and are the primordial 
scalar spectral index and amplitude (at k = 0.05Mpc~^). We marginalize over 
SZ amplitude. The parameters above are given fiat priors, unless otherwise is 
explicitly stated. 

We mostly use two different combinations of cosmological data sets; WMAP5 
data only, and WMAP5 plus other cosmological data (referred to as WMAP5-I--I-) 
For the sake of comparison, in a specific case, we also consider the option of 
omitting BAO constraints, as specified below. 

We then apply priors on v mass according to the KKDG claim and the 
prospected KATRIN results. Data sets and priors used are described in the 
following. The MGMG chains are run on the Titan cluster at Oslo University. 



4.1. Cosmological data 

Firstly, we tested our models using only the fi ve y ear data from the WMAP 
measurements of the GMB radiation (WMAP5) [aESIH- The WMAP5 data 
are analyzed with the Fortran 90 likelihood code provided with the data release. 

In WMAP5-I--I- we then included the galaxy power spectrum from the 2dF 
survey [31,], SNIa data from the SNLS survey [33], and added gaussian priors 



on the Hubble parameter of /i = 0.72 ± 0.08 from the HST key project 



the 



physical baryon density, ujb = 0.022 ± 0.002 [3J, [35!, |36| inferred from the He 
abundance after Big Bang nucleosynthesis, and BAO data from 37 [. 

The BAO analysis in 37[ includes both the SDSS DR7 data and data from 
2dF. Angular distance measurements are compared at z = 0.2 and z = 0.35, 
by using the distance measure d{z) = rs{zd)/ Dv{z). Here rs{zd) is the co- 
moving sound horizon at the baryon drag epoch (see 37,111]) and Dv{z) = 

1/3 

[(1 -I- z)'^ D\z / H {z)\ , where Da is the angular diameter distance. This pro- 
vides information on the late expansion of the Universe, which is important to 
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constrain DE nature. Following [iill^, from the BAO analysis, the inferred 
likelihood £ of a model can be estimated by —2 In £ oc X^-'-C^-'-X, where 



X = 



d(0.2) -0.1905 
d(0.35) - 0.1097 



30124 -17227 
-17227 86977 



(6) 



4.2. KKDC and KATRIN priors 

In order to include KKDC priors on v mass, we use the procedure in refs. [iol 
kill . To account for the dispersion in nuclear matrix {Cmm) estimates, we follow 
[41|, where upper and lower extremes from a compilation of reliable theoretical 
estimates are used, defining a 3cr uncertainty width of Cmm- Combined with the 
uncertainty in T°J^, this results in a KKDC prior of logiQ(TO^^/eV) = —0.23 ± 
0.07 (at Icr). In this mass range the i/ mass eigenvalues are almost degenerate, 
so that we can use that Mi, = 'impp without any loss of accuracy. 

As far as KATRIN is concerned, we use their expected uncertainty (T„2 = 

0.025 eV^ for a Gaussian distribution around a best-fit value m|, as in refs. 
SHIil- We then assume a fiducial value = 0.3 eV (M^ = 0.9 eV), as a 
compromise between cosmological upper limits and KKDC lower limit. 

KKDC and KATRIN priors are imposed in the post processing of MCMC 
chains, using a modified version of the GetDist program provided with the 
CosmoMC package. 



5. Results and discussion 

A first point we want to remark is that including the BAO prior in the 
likelihood analysis causes just minor shifts on the limits, as shown in Table [1] 
and does absolutely not conflict with previous findings. 



Datasets (eV) ^ 

WMAP5 1.70 (1.52) 0.136 ± 0.085 (0.098 ± 0.069) 

WMAP5++ (without BAO) 1.17 (1.13) 0.104 ± 0.047 (0.100 ± 0.043) 

WMAP5++ 1.19 (1.19) 0.094 ± 0.041 (0.092 ± 0.042) 

Table 1 : We report 95% marginalized upper limits on (first column) as well as mean value 
and la uncertainties on /? for the SUGRA (RP) coupled models and different combinations 
of data sets (second column). 

Let US then consider the effects of adding mass priors. 

In Figure [1] we show marginalized 68% and 95% confidence contours in the 
— (3 plane with and without a KKDC prior on Mi,. The left (right) panel 
concerns the SUGRA (RP) potential. Before including KKDC priors, results 
are similar to those in our previous paper Q. In the WMAP5-I--I- case, however, 
where we now include BAO data, the curves shift slightly farther from the 0-0 
case, bringing the 0-0 point well outside Icr contours for both RP and SUGRA 
cases. However, without additional priors on M^, the 0-0 model is not excluded 
by cosmological data alone. 
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Figure 1: 68% and 95% confidence intervals in the — j3 plane. The left panel shows 
resulting limits when using a SUGRA potential, while the right panel shows the limits when 
using a RP potential. Black, thin lines indicate the results when using WMAP5 as the only 
cosmological data set, and red, thick lines show the results when using WMAP5H — h- Dotted 
lines are the cosmology only limits. The resulting limits when also including the KKDC limit 
on are shown with solid lines. We notice that the intersection between /3 and Mi, allowed 
areas, when KKDC controversial results are allowed or omitted, does not vanish but is rather 
small. 



When we add the KKDC prior on M^, we have a significant offset from the 
0-0 modeh ^ = is excluded with 7.3 a (8.0cr)in the SUGRA (RP) case. This 
conclusion relies on the use of WMAP5++; CMB spectra are not sufficiently 
affected by M^, to allow for a statistical detection of /3 alone, even when the 
KKDC prior is imposed. This is not unexpected, as CMB data primarily probes 
the Universe at high z, before DE became important. 

Figure [2] shows ID likelihood distributions for both SUGRA and RP mod- 
els, when using WMAP5-f-+ and the KKDC prior. Here we also plotted the 
mean likelihood for each bin, in addition to the marginalized likelihood. In 
ref. [il a problem that was discussed was the discrepancy between these two 
distributions for the /? parameter, which was caused by a very non-gaussian 
correlation between the log^g A and /? parameters, especially for small values of 
p. From Figure [2] we see that this problem is not very pronounced here. When 
(3 is forced to high values by the KKDC prior, the two probability measures 
correspond quite well to each other, even though we still see some minor shifts. 

In Figure |3] the 2D marginalized likelihood contours are shown when the 
KKDC prior is replaced by the prospected prior from KATRIN with a fiducial 
neutrino mass of = 0.9eV. Although the best-fit here is smaller than 
in the KKDC case, we still get a significant preference for a non-zero /? when 
using WMAP5-h-h. For the SUGRA potential, ^ = gets excluded with 3.9a 
significance, while the corresponding number is 3.6cr for the RP potential. 

The controversial results of the HM experiment, as well as the possibility 
that KATRIN detects a ly mass value, had already triggered a discussion on the 
way to soften cosmological constraints. The best option put forward, up to now, 
was perhaps the possibility of allowing the state parameter of DE, w, to delve 
into the phantom regime. 
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Figure 2: ID likelihood distributions for /3, using WMAP5++ sets and the KKDC prior on 
M,y. Thick, red lines corresponds to the SUGRA model, while thin black lines refers to the the 
RP model. Solid lines denote marginalized likelihood. Dotted lines indicate average likelihood 
of the samples in each bin. 




Figure 3: The same as Figure [T] but using a KATRIN prior with a fiducial neutrino mass of 
= 0.3eV {M^ = 0.9eV) instead of the KKDC prior. Let us also point out that, at variance 
from the KKDC case, the KATRIN prior is fully consistent with cosmological constraints, 
essentially leading to a restriction of the allowed area in the M^-/3 plane. 



Before concluding this Section it is then worth comparing our results with 
those obtainable if we just allow for w < —1, an option already deepened by (l9j . 
Figure |4] shows the likelihood distributions on the My-w plane with or without 
the KKDC constraint. This figure confirms the findings of [13), also when the 
more limited data set considered by them (including, i.e., just WMAP3 outputs) 
is replaced by the whole system of data considered in the rest of this work. 

As is known, SUGRA (or RP) uncoupled cosmologies yield no significant 
likelihood improvement (or worsening) in respect to ACDM. They are therefore 
analogous benchmarks for model likelihood confrontation. It is then clear that, 
if we just allow w to run in the phantom range, we hardly gain a ^ 2-a likelihood 
improvement. On the contrary, if we open the coupling option, already with 
My = 0.9 eV, we approach 4-(t's; passing from 2 to 4-(j means achieving 
statistical significance. 

However, even letting apart statistical evaluations, direct inspection shows 
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Figure 4: Likelihood distribution with or without the KKDC prior (sohd or dotted lines, 
respectively, indicating 1 and 2 cr's), when a DE state parameter u> < — 1 is allowed and 
considering all cosmological data. Notice that the high likelihood areas, with or without the 
KKDC prior, have just a minor intersection at the 2-a level. Let us recall that, with the 
coupling option, there is a (small) intersection even at 1— cr. 



that spectral distortions due to increasing Mj^, both in Ci and P{k), are just 
opposite to those due to (3. On the contrary, w < —1 spectral changes exhibit 
a different scale dependence; in fact, their capacity to allow higher v masses is 
substantially related to their favoring greater 57„ values. 

6. Summary and conclusions 

We have studied the effects of coupling between a dynamical DE component 
and CDM. The observational effects of such a coupling are almost opposite to 
those caused by massive neutrinos, which results in a strong degeneracy between 
the coupling parameter (3 and the neutrino mass M^. 

This suggests the possibility that interactions within the dark sector have 
been hidden to observations, up to now, by the existence of a significant v mass 
value. The v mass, possibly responsible for this chamaleontic effect, lays just 
below the range already inspected in H'^ /3-decay experiments. Such a range 
will be however soon explored by the KATRIN experiment, while the v mass 
detection claimed by a part of the HM team is just above the suitable range. 

In this work, such a possibility was also tested including, for the first time, 
BAO measurements. 

Previous analysis had already shown that a cosmology with mass and 
coupling is statistically preferred to ACDM. The inclusion of BAO data slightly 
increases such preference which, however, still keeps 0(2(7). 

Such /3 — Ml, "degeneracy" can be broken by an external prior on AIi, from 
earth based experiments. If we assume the KKDC claim on neutrino mass to be 
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correct, this then results in a 7-8 a detection of a non-zero /?. If the upcoming 
KATRIN experiments confirm a neutrino mass in the range allowed by the 
KKDC experiment, this wiU on its own standing give a statisticaUy significant 
detection of a non-zero /3. 

Other options considered before, to try to reconcile cosmic data with v mass 
in such range, are far less effective. 

In particular, even if the option w < —1 is added to ACDM models, a 
prior Ml, ~ 1 eV badly modifies the likelihood distribution, so indicating an 
apparent conflict with cosmological measurements. On the contrary, when such 
a prior is added to cosmological data, within the context of /3 ^ models, 
already (slightly) favored by data in respect to ACDM, it just leads to a further 
restriction of the allowed parameter area. In a sense, it appears a "welcome" 
new limit, just narrowing parameter error bars. 

A detection of a non-zero mfj > 0.3eV by KATRIN would radically renew 
the cosmic scenario: the concordance cosmology, ACDM, would be statistically 
falsified; a CDM-DE coupling would become a concrete option. 

One should also keep in mind that in the coming years, with data from the 
Planck satellite, the CMB measurements will improve vastly, further narrowing 
the uncertainties in the (3 — plane, as was shown in 

It may be then worth pointing out that admitting a linear interaction be- 
tween CDM and DE, as done here, is just the next approximation to assuming 
them fully decoupled. The idea underlying this analysis is that modeling the 
Dark Cosmic Sector with two independent components is just a first step to- 
wards understanding its complex nature. After a hypothetical confirmation of 
the coupling option by laboratory data, further cosmic data might allow to go 
even beyond the assumption of linear interaction, possibly providing a more 
detailed insight into one of basic question of modern physics, the nature of the 
cosmic dark components. 
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